APERIODIC FRACTIONAL OBSTACLE PROBLEMS 



MATTEO FOCARDI 

Abstract. We determine the asymptotic behaviour of (bilateral) obstacle problems for fractional 
energies in rather general aperiodic settings via F-convergence arguments. As further developments 
we consider obstacles with random sizes and shapes located on points of standard lattices, and the 
case of random homothetics obstacles centered on random Delone sets of points. 

Obstacle problems for non-local energies occur in several physical phenomonona, for which our 
results provide a description of the first order asympotitc behaviour. 



1. Introduction 

Non-local energies and operators have been actively investigated over recent years. They arise in 
problems from different fields, the most celebrated being Signorini's problem in contact mechanics: 
finding the equilibria of an elastic body partially laying on a surface and acted upon part of its 
boundary by unilateral shear forces (see [43], [3f]). In the anti-plane setting the elastic energy can be 
then expressed in terms of the seminorm of a H^/^ function, or equivalently as the boundary trace 
energy of a W^^'^ displacement. 

As further examples we mention applications in elasticity, for instance in phase field theories for 
dislocations (see [.34] and the references therein); in heat transfer for optimal control of temperature 
across a surface [33], [G]; in equilibrium statistical mechanics to model free energies of Ising spin 
systems with Kac potentials on lattices (see [2] and the references therein); in fluid dynamics to 
describe flows through semi-permeable membranes [30]; in financial mathematics in pricing models 
for American options [4]; and in probability in the theory of Markov processes (see [10], [11] and the 
references therein). 

Many efforts have been done to extend the existing theories for (fully non-linear) second order elliptic 
equations to non-local equations. Regularity has been developed for integro-differential operators 
[10], [11], [22], and for obstacle problems for the fractional laplacian (see [20], [44] and the references 
therein). Periodic homogenization has been studied for a quite general class of non- linear, non-local 
uniformly elliptic equations [41] and for obstacle problems for the fractional laplacian [19], [32]. 
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In this paper we determine the homogenization hmit for bilateral obstacle problems involving frac- 
tional type energies. In doing that we employ a variational approach by following De Giorgi's F- 
convergence theory. Many contributions in literature are related to the analogous problem for ener- 
gies defined in ordinary Sobolev spaces or equivalently for local elliptic operators. The setting just 
mentioned will be referred to in the sequel as the local case in contrast to the non-local framework 
object of our analysis. Starting from the seminal papers by Marchenko and Khruslov [ Ranch and 
Taylor [39] , [40] , and Cioranescu and Murat [24] there has been an outgrowing interest on this kind of 
problems with different approaches (see the books [7], [15], [IG], [23], [28] and the references therein 
on this subject). We limit ourselves to stress that F-convergence theory was successfully applied to 
tackle the problem and to solve it in great generality (see [29], [26], [27]). 

Let us briefly resume the contents of this paper in a model case (for all the details and the precise 
assumptions see section 3). With fixed a bounded set T, and a discrete and homogeneous distribution 
of points A = {x^jigz" (see Definition 2.2), define for all j e N the obstacle set Tj C R" by 
Tj = UigZ" \£j 4- e"^^" *p) 2^ j ^ -^j^cj-q (e^ )jg]vj is a positive infinitesimal sequence. Given a bounded. 



open and connected subset U of R" , n > 2, with Lipschitz regular boundary we consider the functionals 
J'j : L'l'iU) [0, +oo] given by 



and -l-cxD otherwise. Here, W^''p{U) is the Sobolev-Slobodeckij space for s S (0,1), p <£ (l,4-oo) and 
sp e (1, n), capj p is the related variational (p, s)-capacity, and u denotes the precise representative of 
u £ W^'P{U) which is defined except on a cap^ p-neghgible set (see subsections 2.4 and 2.5). 

In Theorem 3.3 we show that the asymptotic behaviour of the sequence (.?^j")jeN is described in 
terms of F(LP(J7))-convergencc by the functional JF : L'p{U) [0, 4-oo] defined by 



if w e W^^'P([/), 4-00 otherwise in LP{U). The quantities 6 and /3 represent respectively the limit 
density and the limit normalized distribution of the set of points SjA and can be explicitely calculated 
in some cases (see (3.1) and (3.3) for the exact definitions. Remark 3.2 for further comments, and 
Examples 3.5-3.7 where some cases are discussed in details). 

Adding zero boundary conditions, F-convergence then implies the convergence for minimizers and 
minimum values of (jF,)jgN to the respective quantities of J^. More generally, we can study the 
asymptotic behaviour of anisotropic variations of the Gagliardo seminorm (see subsection 3.4). 

We remark that it is not our aim to establish a general abstract theory as Dal Maso did in the local 
case [26], [27] (see [!)] for related results in random settings), nor to consider the most general framework 
for homogenization as the one proposed by Nguetseng in the linear, local and non degenerate case 
[38] ; but rather we aim at giving explicit constructive results for a sufficiently broad class of fractional 
energies and non-periodic obstacles. 




r \u{x)-u{y)\P 
uxu \x-y\"+'P 



dxdy if u e ^'^{11), u = Q cap^ ^ q.e. on Tj n U 
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The main novelty of the paper is that we extend the asymptotic analysis of obstacle probelms for local 
energies to non-local ones. In doing that we use intrinsic arguments and give a self-contained proof. 
We avoid extension techniques with which the original problem is transformed into a homogenization 
problem at the boundary by rewriting fractional seminorms as trace energies for local (but degenerate) 
functionals in one dimension higher (see [25], [19], [32]). This is the well known harmonic extension 
procedure for H^/'^ functions; more recently it was proved to hold true also for the fractional laplacian, 
corresponding to p = 2 and s G (0, 1) above, by Caffarelli and Silvestre [21]. Because of our approach, 
we are able to deal with non-local energies which can not be included in the previous frameworks. 
Generalizations are likely to be done in several directions, the most immediate being the choice of the 
obstacle condition. For, we have confined ourselves to the basic bilateral zero obstacle condition in 
the scalar case only for the sake of simplicity, improvements to vector valued problems and unilateral 
obstacles seem to be at hand (see [5], [32]). 

Our analysis do not need the usual periodicity or almost periodicity assumptions for the distribution 
of obstacles. We deal with aperiodic settings defined after Delone set of points (see subsection 2.3), 
the set A introduced above. No regularity or simmctry conditions are imposed on A, only two simple 
geometric properties are assumed: discreteness and homogeneity. They turn out to be physically rea- 
sonable conditions; as a matter of fact Delone sets have been introduced in ?i-dimensional mathematical 
cristallography to model many non-periodic structures such as quasicrystals (see [42]). Essentially, 
these assumptions guarantee that points in A can neither cluster nor be scattered away. 

Furthermore, we show that with minor changes the same tools are suited also to deal with some 
random settings. In particular, we deal with obstacles with random sizes and shapes located on points 
of a standard lattice in R", a setting introduced by Caffarelli and Mellet [18], [19]; and consider also 
homothetics random copies of a given obstacle located on random lattices following Blanc, Lc Bris 
and Lions [13], [14] (for more details see section 4). 

As a byproduct, our approach yields also an intrinsic proof of the results first obtained by [19] 
that avoids the extension techniques in [2]]. A different proof using F-convergence methods, but still 
relying on those extension techniques, was given by the author in [32]. 

The key tools of our analysis are Lemmas 3.8 and 3.9 below. By means of these results we reduce 
the F-liniit process to families of functions which are constants on suitable annuli surrounding the 
obstacle sets. Lemma 3.9 is the counterpart in the current non-local framework of the joining lemma 
in varying domains for gradient energies on standard Sobolcv spaces proven by Ansini and Braides 
[■^)]. It is a variant of an idea by De Giorgi in the setting of varying domains, on the way of matching 
boundary conditions by increasing the energy only up to a small error. As in the local case the proofs 
of Lemmas 3.8 and 3.9 exploit De Giorgi's slicing/averaging principle and the fact that Poincare- 
Wirtinger inequalities are qualitatively invariant under families of biLipschitz mappings with cqui- 
bounded Lipschitz constants. 
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Despite this, the non-local behaviour of fractional energies introduces several additional difficulties 
into the problem: Lemmas 3.8 and 3.9 do not follow from routine modifications of the arguments used 
in the local case. New ideas have to be worked out mainly to control the long-range interaction terms. 
A major role in doing that is played by the counting arguments in Proposition 2.5, Hardy inequality 
(see Theorem 2.8), and the estimates on singular kernels in Lemma A.l. 

The paper is organized as follows. In section 2 we list the necessary prerequisites on F-convergence, 
Sobolev-Slobodeckij spaces and variational capacities giving precise references for those not proved. 
Section 3 is devoted to the exact statement and the proof of the homogenization result for deterministic 
distribution of obstacles. To avoid unnecessary generality we deal with the model case of fractional 
seminorms. Generalizations to anisotropic kernels are postponed to subsection 3.4. The ideas of 
section 3 are then used in section 4 to deal with the two different random settings mentioned before. 
Finally, we give the proof of an elementary technical result, though instrumental for us, in Appendix A. 

2. Preliminaries and Notations 

2.1. Basic Notations. We use standard notations for Lebesgue and Hausdorff measures, and for 
Lebesgue and Sobolev function spaces. 

The Euclidean norm in R" is denoted by | • |, the maximum one by | • |oo. Br{x) stands for the 
Euclidean ball in R" with centre x and radius ?' > 0, and we write simply in case a; = 0. As usual 

Given a set E C R" its complement will be indifferently denoted by E"^' or R" \ E. Its interior 
and closure are denoted by int(E) and E, respectively. Given two sets E CC F in R", a cut-off 
function between E and F is any € Lip(R", [0, 1]) such that ip\-^ = 1, '^\-R,n\F = 0, and Lip((^) < 
l/dist(£;,9F). 

Given an open set A C R" the collections of its open, Borel subsets are denoted by A{A), B{A), 
respectively. The diagonal set in R" x R" is denoted by A, and for every (S > its open (5-neighborhood 
by := {{x,y) e R" x R" : \x — y\ < 5}. Accordingly, for any set E C R" and for any J > 

Es e R" : dist(x, E) < 8}, E^s := {x e E : dist(a;, dE) > S}. (2.1) 

In the following U will always be an open and connected subset of R" whose boundary is Lipschitz 
regular. 

In several computations below the letter c generically denotes a positive constant. We assume 
this convention since it is not essential to distinguish from one specific constant to another, leaving 
understood that the constant may change from line to line. The parameters on which each constant 
c depends will be explicitely highlighted. 
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2.2. F-convergence. We recall the notion of r-convcrgcncc introduced by De Giorgi in a generic 
metric space {X,d) endowed with the topology induced by d (see [28], [15]). A sequence of func- 
tionals Fj : X [0, +oo] T-converges to a functional F : X ^ [0, +oo] in u G X, in short 
F{u) = r-linij Fj{u), if the following two conditions hold: 

(i) (liminf inequality) V (uj)jeN converging to u in X, we have liminfj Fj{uj) > F{u); 

(ii) {limsup inequality) 3 (uj)jgN converging to u in X such that limsupj Fj{uj) < F(u). 

We say that Fj T-converges to F (or F= T-limjFj) if F{u) = T-limj Fj{u) Vu G X. We may also 
define the lower and upper T-limits as 

F- lim sup i^j (u) = inf{limsup-Fj(uj) : Uj — u}, 

3 3 

V- lim inf Fj (u) = inf {lim inf Fj (uj) : Uj —^ u}, 

3 3 

respectively, so that conditions (i) and (ii) are equivalent to F-limsup ji^j(u) = F-liminfji^j(u) = F{u). 
Moreover, the functions F-limsup_,-Fj and F- liminf ji^j are lower semicontinuous. 

One of the main reasons for the introduction of this notion is explained by the following fundamental 
theorem (see [28, Theorem 7.8]). 

Theorem 2.1. Let F = F-linij Fj , and assume there exists a compact set K <Z X such that mix Fj = 
vaiK Fj for all j . Then there exists minx = limjiidxFj. Moreover, i/ (uj)jgN is a converging 
sequence such that limj Fj{uj) = linij infx Fj then its limit is a minimum point for F. 

2.3. Non-periodic tilings. In the ensuing sections we will deal with a general framework extending 
the usual periodic setting. The partition of R" we consider is obtained via the Voronoi' tessellation 
related to a fixed Delone set of points A. We refer to the by now classical book of M. Senechal [42] 
for all the relevant results. 

Definition 2.2. A point set A C R" is a Delone (or Delaunay) set if it satisfies 

(i) Discreteness.' there exists r > such that for all x,y G K, x ^ y, \x — y\ > 2r; 
(ii) Homogeneity or Relative Density; there exists R > such that An Bfi{x) ^ for all x £ R". 

It is then easy to show that A is countably infinite. Hence, from now on we use the notation 
A = {x^jigz"- By the very definition the quantities 

rA:=iinf{|a--y|: x,j/eA, x^?/}, i?A := mf{i? > : A n Bi^(x) ^ Vx G R"} (2.2) 

are finite and strictly positive; i?A is called the covering radius of A. 

Definition 2.3. Let A C R" he a Delone set, the Voronoi cell of a point G A is the set of points 

V' := {y e R" : |y - x^l < 12/ - x'^j, for all i ^ k}. 
The Voronoi tessellation induced by A is the partition o/R" given by {V^^jigz"- 
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In the following proposition wc collect several interesting properties of Voronoi tessellations (see [42, 
Propositions 2.7, 5.2]). 

Proposition 2.4. Let A C R" be a Delone set and {T^^jiez" Us induced Voronoi tessellation, then 
(i) the V^'s are convex polytopes fitting together along whole faces, and have no interior points 
in common; 

(ii) if and share a vertex z, then and lie on dBp{z) with An Bp{z) = 0, p < i?A- 
Hence, {F^jigZ" is a tiling, i.e. the 's are closed, have no interior points in common and UiV^ = 
R". More precisely, 

(Hi) {V^}ii=z" is a normal tiling: for each tile we have Brf^{x^) C V^^ C i?/f^(x^); 
(iv) {V^}i^z" is a locally finite tiling: #(A n Bp{x)) < +oo for all x € R", p > 0. 

Further properties that will be repeatedly used in our analysis are summarized below. We omit 
their proofs since they are justified by elementary counting arguments. For any A G ^(R") we set 

Xa(A) := {i e Z" : t/^ C A}, ^a(A) {i e Z" : n 9/1 7^ 0}. (2.3) 

Proposition 2.5. Let A C R" be a Delone set and {y^}igz" its induced Voronoi' tessellation. Then, 

cOrymMA)) < £"(A), u;„rU{.yA{A)) < {dA)n,, £" {A\Uj,^a)V') < {dA)n,. (2.4) 
In particular, there exists a constant c — c{n) > such that for every i £ Z", m € N if holds 

#{k e Ia{A) : mrA < - x'^U < {m + l)rA} < cm''-\ (2.5) 

2.4. Sobolev-Slobodeckij spaces. Let A C R" be any bounded open Lipschitz set, p e (l,+oo), 
s e (0, 1) and ps e (!,«), by ^'^{A) we denote the usual Sobolev-Slobodeckij space, or Besov space 
Bpp{A). The space is Banach if equipped with the norm ||m||m/s.p(a) — \\u\\lp{A) + |^t|vy=.j'(A); where 

p f \u{x) - u{y)\P 

■■=J,^, \,.y\n,s, d^dy- 

We will use several properties of fractional Sobolev spaces, giving precise references for those employed 
in the sequel in the respective places mainly referring to [1] and [45]. 

In the indicated ranges for the parameters p, s it turns out that W'^ is a reflexive space (see [45, 
Thm 4.8.2]), Sobolev embedding and Sobolev-Gagliardo-Nirenberg inequality hold (see [1, Chapter 
V]), and traces are well defined (see (2.9) below). We remark that the restriction on s are necessary, 
since otherwise W'''P{A) contains only constant functions if s > 1, while for ps < 1 traces are not well 
defined (see also (2.9) below). The exclusion of the other cases is related to capacitary issues (see 
subsection 2.5). 

Poincare-Wirtinger and Poincare inequalities in fractional Sobolev spaces are instrumental tools in 
the sequel. Thus, we state explicitely those results in the form wc need. Their proof clearly follows 
from the usual argument by contradiction once W^'^ is reflexive and is endowed with a trace operator. 
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Theorem 2.6. Let n > 1, p Cz (l,+oo), and s e (0,1). Let A C R" be a bounded, eonnected open 
set, and O any measurable subset of A with C"{0) > 0. Then for any function u e W^''p{A), 

11" - '"oIIlp(^) < cpw|w|^.,p(^), (2.6) 
for a constant cpw = cpw{n,p, s,0, A). Moreover, for any u G Wq''^{U) we have 

||u|liP(A) < cp|u|^„p(^), (2.7) 

for a constant cp — cp{n,p, s, A). 

Remark 2.7. Let {^t)t£,3r be a family of biLipschitz maps on A with sup5.(Lip($t) + Lip{^^^)) < 
+00, then a simple change of variables implies that the constants cpwin,p, s,^tiO),^tiA)) and 
cp{n,p, s, $t(A)) are uniformly-bounded in t. In particular, a scaling argument and Holder inequality 
yield for any z G R" and r > and for some c = c{n,p, s, O, A) > 

\\u - u,+ro\\l,(^,+,A) ^ cr^P|"l^^=.P(.+.A)- (2-8) 
A similar conclusion holds for Poincare inequality (2.7). 

Next, we recall the fractional version of Hardy inequality (see [45, Theorem 4.3.2/1, Remark 2 pp. 
319-320] and [40] for further comments). To this aim we introduce the space W^'''^{Bi) := {u E 
W'-P(R") : sptu C Bi}. It is clear that C^{Bi) is dense in W''P{Bi) for any p e (1, +oo) and that 
W'''P{Bi) C Wq'''{Bi). The latter inclusion is strict if s — 1/p G N, and more precisely it holds 

{W'^PiBi) if s > 1/p- 1, s - 1/p^ N 
W^PiBi) if s e (0, 1/p]. 

Theorem 2.8. There exists a constant c = c{n,p, s) such that for every u G W'^'P{Bi) we have 

\u{x)\P 



Ib^ dist(a;,aSi)«P 

Remark 2.9. The usual scaling argument, Poincare inequality (2.7) and Theorem 2.8 then yield 

\u{x)\P 



/^„dist(x,as.-F^^-^"""'^-^(^^) 

for every r > and u G W'''P{Br), for a constant c = c{n,p, s) independent from r. 

2.5. Fractional capacities. We recall the notion of variational capacity for fractional Sobolev spaces 
and prove some properties relevant in the developments below. Those properties, straightforward in 
the local case, require more work in the non-local one. Let p G (1, +oo) and s G (0, 1) be given as 
before, for any set T C R" define 

cap,„(r):= inf inf ( lul^, p,„„, : u G iy'^P(R"), u > 1 £" a.c. on , (2.11) 
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with the usual convention inf = +00. The set function in (2.11) turns out to be a Choquet capacity 
(see [1, Chapter V]). Recall that a property holds cap^ ^ quasi everywhere, in short capj, ^ q.e. on A, if 
it holds up to a set of cap^ ^ zero. In particular, any function u in W'''P{A), A £ ^(R"), has a precise 
representative u defined cap^ ^ q.e. and the following formula holds (see [1, Proposition 5.3]) 

cap,_^(T) inf : w G M^^'^'(R"), w > 1 q.e. on t} . (2.12) 

Coercivity of the fractional norm is ensured only in the (R") topology, p* :~ np/{n — sp) is the 
Sobolev exponent relative to p and s (see [1, Chapter V]). Thus, a minimizer for the capacitary problem 
exists in the homogeneous space M^*'P(R") = {u E (R") : |u|vi/s,p(-R^n-) < +00}. Uniqueness is 
guaranteed by the strict convexity of the fractional energy, thus the minimizer of (2.11), (2.12) will 
be denoted by and called the capacitary potential for T. 

Remark 2.10. If ps > n points have positive capacity and W^'^ is embedded into C*^. In this case 
it is well known that the homogenized obstacle problem turns out to be trivial, this is the reason why 
we required ps < n. The borderline case ps = n deserves an analysis similar to that we will perform 
but different in some details, so that its study is not dealt with in this paper (see for instance [2-1] and 
[17] in the local framework). 

Let us prove that set inclusion induces an ordering among capacitary potentials. As a byproduct 
we also show that admissible functions in the capacitary problem can be taken with values in [0, 1]. 

Lemma 2.11. IfT C F, then < < < 1 £" a.e. on R". 

Proof. First, take note that for all u, v E ^[^(.(R") we have 

\{u V v){x) ~{uV v){y)\ < \u{x) - u{y)\ V \v{x) - v{y)\, 
\{u A v){x) - (u A v){y)\ < \u{x) - u{y)\ V \v{x) - v{y)\. 

In particular, uniqueness of the capacitary potential U"^ and by choosing w = in (2.13)i and v = 1 
in (2.13)2 imply that < < 1 £" a.e. on R" for any subset T C R". 

Moreover, (2.13) yields that V and A e iy*'P(R"). Set Ut = {u^ < u^} and 
Up = {u^ < u'^}, and assume by contradiction that C"'{Uf) > 0. By taking V as test function 
in the minimum problem capg p{F) and recalling the strict minimality of , we infer Im^I^s.p^r^™) < 
\u'^ V u^l^s.p^j^n-) . An easy computation then leads to 

+ 2 / ^"^y-^y ^xdy < n.,^,,^, + 2 / ^""y-fjf^ dxdy. 

^ JutxUi. \x-y\''+'P JutxUj. |a:-yr+"P 

The latter inequality can be used to estimate the fractional norm of u'^ A as follows 



1"^ A ^^^|^.,p,R„^ < k^|L..p.H„^ (2.14) 



r \u^{x) ~ u''iy)\P + lu^jx) ~ u^{y)\P - \u^{x) - u^{y)\P ~ \u^{x) ~ uP{y)\P 



dxdy. 
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We claim that the second term on the rhs in (2.14) is non-positive, this would imply that the strict 
inequality sign holds above, which in turn would give a contradiction since v?" A is a test function 
for the capacitary problem related to T. 

To conclude consider the auxiliary function '0(s, t) := \s — u'^{y)\P + \u'^{x) — t\P —\s — t\P — lu^iy) — 
u'^{x)\P and the set H := {{s,t) : s > u^{x), t < u^{y)}. Elementary calculations yield that 
maxH ^ < 0- Take note that the numerator of the integrand of the second term in the rhs of (2.14) 
equals to ip{u^{x),u^{y)) and u^{x) > u^{x) £" a.e. x G Ut, and u^{y) < u^iy) -C" a.e. y e Up. □ 

In the sequel we arc interested into relative capacities, for which we introduce two different notions. 
The first one is useful in the F-liminf inequality, the second in the F-limsup inequality, respectively. 
For every < r < R set 

cap,^p(T,Bfl;r) inf : w S W''P(R''), i« = on R" \ B^, w > 1 q.e. on t} , 

and 

Cs,p{T,Br) := m{^\w\p^^,,^^„^ : w G M^^'^'(R"), w = on R" \ Si?, w > 1 q.e. on t} . 

To prove the convergence of relative capacities to the global one we introduce some notation to simplify 

the calculations below: for any £"-measurable function w and any "-measurable set E <Z U x U 

consider the locality defect of the W^'P scminorm 

^ / f \w{x) — w{y)\P , , 

VsA^,E) :^JJ-j-L-^dxdy. 

The terminology is justified since given two disjoint subdomains A, B (Z U, it holds 

\MW'.P{AUB) = \^\w^.piA) + \Mw^.P{B) + '^-T^sA'^, A X B). (2.15) 

In particular, Vs^piw^A x A) = \w\w''P{A)- 

We are now in a position to prove the claimed converge result for relative capacities. Actually, we 
show uniform convergence for families of equi-bounded sets. A different argument yielding pointwise 
convergence will be exploited in the generalizations of subsection 3.4 (see Lemma 3.12). The latter 
is sufficient for the proof of Theorem 3.3; the advantage of the approach below is that it can be 
carried over straightforward to the case of obstacles with random sizes and shapes for which uniform 
convergence is necessary (see subsection 4.1). 

Lemma 2.12. For all p > it holds 

lim sup |C,,p(r,B,)-cap_(r)| =0. (2.16) 

Moreover, there exists a constant c = c(n, s,p) such that for all < p < r < R 

cr^P 

sup (cap,.p(T) - cap,.p(T, Br- r)) < — —C^ABp, Br). (2.17) 

TCB„ — ryp 
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In addition, if R{r)/r +00 as r +00 we have 

lim sup |cap,^p(T)-cap,^p(T,Sfl,(,);r)| =0. (2.18) 

Proof. It is clear from the very definitions that (0, +00) 3 r ^ Cs,p{T,Br) is monotone decreasing, 
and moreover that cap, ^(T, Bb,\ r) < Cs.p{T^ Br)-, cap, p(T) < Cs.p{T, Br). 

Let us first prove (2.16). To this aim take , u^p £ M^'''^(R") the capaeitary potentials of T 
and Bp, respectively; by Lemma 2.11 we have < m"'" < u^^ < 1 £" a.e. on R". In addition, u^p 
is radially symmetric and decreasing (to zero) since the fractional norm is stricly convex, rotation 
invariant and decreasing under radial rearrangements (for the last result see [A, Section 9]). 

With fixed 6 > consider the Lipschitz map ipsit) ■= jE^ V 0, hip^ips) < (1 — ^^nd set 

ws{x) :— il}s{u^{x)). Up to £" negligible sets, {ws > 0} = {u^ > 5} Q {u^" > 5] Q Br^, for some 
Rs +00 as (5 ^ 0+. Then ws G VK"'?'(R") with 

Moreover, wg > 1 q.e. on T. Since cap, „(•) is an increasing set function, we infer 



< CsAT.Rs) - cap,^p(r) < (^(Y^ " ^ ) ^^P«^p(^) ^ KU^sy " ^ ) '^p^.p^^'')- 

In conclusion, (2.16) follows from the monotonicity properties of r ^ Cs,p{T,Br). 

To prove (2.17) take note that any admissible function u for the minimum problem defining eap^ ^(r, Br 
is admissible for the one defining cap^ piT), too. Furthermore, for some constant c = c(n, s,p) it holds 



Br 



, ,,, c 

< 



In the last inequality we used the scaled version of Poineare inequality (2.7) as follows from Remark 2.7. 
By passing to the infimum on the admissible test functions we infer 

cap,jT)-cap,jT,Bn;r) < ^^^—^cap^^^iT, BR-r) < ^__-^C.,p(T, i?,). 

We deduce statement (2.17) since Cs,p{-,Br) is a monotone increasing set function. 
Eventually, (2.18) follows at once from (2.16), (2.17), and the fact that cap, ,p(T, Br; r) < Cs,p{T, Br). 

□ 



Remark 2.13. Clearly estimate (2.17) blows up for r ^ R. In such a case by using Hardy inequality 
one can only prove that capg p{T) < (1 + c)cap, ^(r, Br, R) for some c = c{n, s,p) > 0. 
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Remark 2.14. If S,r is a {l/r)-minimizer for Cs,p{T,Br) and R{r)/r +00 as r +00, then 

lim T>s.p{^r,Bii(^r) X = 0. 

Indeed, being admissible for the problem defining cap^ p(r, Bj^; r), for all r < R, with 

from (2.18) we conclude. 

3. Deterministic Setting 

3.1. Statement of the Main Result. Consider Delone sets Aj = {x^jigz", and let Rj := Ra^, 
Ij{A) := lA.iA), yj{A) := J?a^(A), for all A e AiU) (see (2.3) for the definition of I,, y^). Fix 
rj S (0, TAj], and assume that the r^'s and A^'s are such that 

limrj = 0, (1 <) limsup(i?j/rj) < +00, (3.1) 

3 j 

lim#X,(C/) r" = e e (0, +00), (3.2) 

'"■'■^=^x7Zn ^ -^xj L[/ u;*-a(C/), (3.3) 
for some /? G i^([/, [0,+oo]) with ||/3||l1(c/) = 1- 

Remark 3.1. Condition (3.2) implies that rj ^ r^j since limsupj #X, ([/) rj^^ < +00 by (2.4)i and 
liminfj ij=I.j{U)r\^ > by (2.4)2 ("see aZso Remark 3.4). 

Remark 3.2. It is well known that the w*-Cb{U) convergence of to /i in (3.3) can be restated 

as 

fij{A) n{A) for all A e A{U) with fi{U n dA) 0. (3.4) 

Proposition 2.5 and conditions (3.1), (3.2) imply that assumption (3.3) is always satisfied up to sub- 
sequences. First, let us show that any w*-Cq{U) cluster point of the probability measures (/ij)jgN is 
absolutely continuous w.r.to C"\—U. For, let fij converge to fj. w* -Cq{U), then since Brj(x^) Q Vj' , 
for every ip £ Cq{U) and S > uniform continuity yields for j sufficiently big 

'^^'^J' = //-^^r^^ ^'^'^j) - jiT^nn 1 fdx + S < }^ [ \Lp\ dx + S. 

U ^^^^^^ I^) "^^3^^' X~iU)''^f (^ryf#Ij{U) Jjj 

By taking as test functions ±ip and first letting j — > +00 and then S ^ 0^ we infer 

\(p\ dx. 

The latter inequality implies fi <C C"\—U. Actually, since Vj' C Bji.(yij) arguing as above it follows 
P e L°°{U) with (liminfj rj/i?^)" < ujn0(3{x) < 1 for a.e. xeU. 

Furthermore, take A' , A G A(U) such that A' CC A CC U , then for j sufficiently big we have 









<—[ 


Ju 


^nO Ju 
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_ #{iGX,(t/): x^€^} 4I^{A') + #J,{A') 



^ _ #UU\A') (2.4),^(2.4), ^ C-{U\A') 



Hence, equi-tightness o/(/^j)jgN follows from 

i„..f,,w>i-,i,„»,,(^^)"£^. 

Thus, Prokhorov theorem gives the w*-Cb{U) compactness in (3.3) up to subsequences. 

With fixed a bounded set T, for all j G N define the obstacle set Tj C R" by Tj = Diez^Tj where 

xj + A,r, and A, := rj/^""^^' . (3.5) 

Take note that Tf C for all i e Z" and j e N. 
Consider the functionals Tj : LP{U) [0, +oo] defined by 

^ , , f if " e W^'PiU), u = cap,.p q.e. on T, n C/ 

•^iW = < ^ (3-6) 

+00 otherwise. 

Theorem 3.3. Let U G ^(R") be bounded and connected with Lipschitz regular boundary; and assume 
that (3.1)-(3.3) are satisfied. 

The sequence (J-'j)jgN T-converges in the V{U) topology to J- : LPiJJ^ — > [0, +cx)] defined by 

:F{u) = +0cap,,p(r)^|«(rr)|P/?(a;)(ix (3.7) 

ifu e VF"'P([/),, +00 otherwise in Lp{U). 

Remark 3.4. If 6 G {0, +00} simple comparison arguments show that the conclusions of Theorem 3.3 
still hold true, though the T-limit is trivial in both cases. 

Let us show some examples of sets of points included in the framework above. In the sequel [ej)j^^<^ 
will always denote a positive infinitesimal sequence. 

Example 3.5. Given a Delone set of points A in R" let Kj :~ SjA, then r^j ~ ^j^A ^.i^-d Raj — £jRA- 
If rj ^ assumptions (3.2) and (3.3) hold true up to the extraction of a subsequence according to 
Remark 3.2, respectively. Several ways of generating Delone sets of points are discussed in [42]. 

More explicit examples can be obtained as follows (see Examples 4.9, 4.10 for the stochastic versions). 

Example 3.6. Let 4> : R" R" be a diffeomorphism satisfying 

||V$||l=o(r.„) < M, and inf det V$ > > 0. 



APERIODIC FRACTIONAL OBSTACLE PROBLEMS 



13 



Then the smallest eigenvalue o/V*<f>V$ is greater than vM^ and thus for all x,y G R' 



n 



iyM^-"\x -y\< \^{x) - <P{y)\ < M\x - y\. 



Set Aj ~ {$(eji)}igz"7 then (vM^ ^ /2)£j < r\^ < < Msj. An easy computation and (3.4) 
yield the w* -Cb{U) convergence of the measures /ij in (3.3) to fi — j3C'^\_U with 



\JU / 
Eventually, if rj e (0, r\.] with rj/sj —^J>Owe have = Ju det V<i>^"'^(a;) dx. 

Example 3.7. Let ^ be a diffeomorphism as in the previous example, and set Aj ~ {ej^{±)}±:^z" ■ 
As before, we have {vM^^'^ /2)£j < r\^ < < Mej, by using (3.4) it can be checked that if 

(det V<i>~"'^)(-/£j)jgN converge to g weakly * in Lj^^(R"), the measures in (3.3) converge w*-Cb{U) 
to fi = I3C^\-U with 



By choosing rj € (0, r^^] with rj/sj j > 0, we have = Jjj g{x) dx. 

3.2. Technical Lemmas. To prove Theorem 3.3 we establish two technical results which are instru- 
mental for our strategy. Roughly speaking we show that the F-limit can be computed on sequences of 
functions matching the values of their limit on suitable annuli surrounding the obstacle sets. To give 
a proof as much clear as possible we first work in an unsealed setting in Lemma 3.8, and then turn 
to the framework of interest in Lemma 3.9. The method of proof is elementary and based on a clever 
slicing and averaging argument, looking for those zones where the energy does not concentrate. The 
relevant property we prove is that the energetic error of the construction wc perform is estimated by 
a local term: a measure. 

Lemma 3.8. Let A be a Delone set of points. For any m G N, m > 2, p G (0,rA/2) and i G X\{U) 
let A[ = + Bp/jn \ Bp/jn'^ , Ai ~ x.^ + Bp \ B p/j^^, and ipi{-) = (y5(- — x^), where Lp is a cut-off function 
between Bp/^a \ Bp/^^i and Bp \ Bp/^3 . 

Then there exists a constant c = c{n,p,s) > such that for any function u G W'^'^{U), and any 
ij^If^{U) -tuple of vectors {zi}^^XK{U)i ^ R"; the function 



belongs to ^'^{U), u = z± on A[ and w = u on U \ A, with A := Uigi^([/)Ai; in addition for every 
measurable set E Q U x U it holds 



\Vs,piw,E)~Vs,piu,E)\<c \v,^piu,UxA)+m^Pp-'P ^ / \u{y) ~ z.f dy \ . (3.8) 
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Proof. By construction w = u on U \ A and u ~ Zi on A'^ for each i. 
To prove that w G W'^'P{U) and estimate (3.8) in case E ^ U x U we use (2.15) to get 

H^w^.PiU) = \<w^..iu\A) + + 2I?.,p(u', AxiU\ A)). (3.9) 

In order to control the last two terms on the rhs above we will use two different splitting of the 
oscillation wlx) — w{y) corresponding roughly to short range and long range interactions estimates. 
First decompose further the seminorm of w on A in (3.9) as follows, 

\^\w'..iA)^ E / ■■■dxdy+ I .■■dxdy=:h+l2. (3.10) 

lelAiU) {(i,k): i^itk} 

Next we deal with Ii: since lySi = on Ai for 1 7^ i and < ip± < 1 we have 

J ^ f 1(1 - 'fi{x))u{x) - (1 - 'fi{y))u{y) + {(pijx) - (pi{y})zif ^ ^ 

ieiA(;7) ^ ./AixAi I y\ 

Take note that Lip((^i) < 2w? j p for all i € Ia{U), then Fubini theorem and (i) in Lemma A.l applied 
with 1/ = n + {s — l)p and O = Ai imply 

- J ji*..... )7.,.., 

<c J2 fH'vK=.nA.)+"^''^''' / \uiy)-z,\Pdy). (3.11) 
To estimate /2 we rewrite w as follows 

w{x) - w{y) = u{x) - u{y) + ^ (pi(x)(zi - w(a;)) + ^ (^i(y)(u(y) - Zi). (3.12) 

1&Ia{U} 1GIa(C/) 

With the help of (3.12) and since tpi = on A \ Ai we infer 

J ^ Sr^ f \u{x) - u{y) + ipi{x){zi - u(a;)) + ipy,{y){u{y) ~ Zy,)\P 

Thus, we can bound each summand in I2 as follows 



\w{x) -w{y) \^ ^^^^ ^ 



dxdy < iP-^Vs,p{u, Ai X A) (3.13) 



.3.-1 / dx / i^iWK^)-^i)r ^^^3P-i / / i^k(^)(u(,)-z.)r ^^^ 



Further, take note that for every fixed (i,k), with i 7^ k, if .t G and y G we have |a; — y| > 
dist(Ai, ^ij) > |a;^ — x^\ — r^, and then \x — y\ > \x^ — x^\/2 since |a;^ — x^\ > 2r\. Hence, being 
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< uJnP^ for all k e Z", from (2.5) and the choice p £ (0, rA/2) it follows for some constant 
c = c{n) > 

/i>l {k^l, hrA<bi-xM=o<(''+l)'-A} '»>i 

By summing up on (i,k), i 7^ k, Fubini theorem, (3.13) and (3.14) imply for some c = c{n,p,s) > 



/2<3J'-il?,.p(w,Ui(^iXU^iA)) + 4r E / |M(a;)-Zi|Pdx, 



(3.15) 

ieiA(;7) ■' ' 

We turn now to the locality defect term in (3.9). We use again equality (3.12) and notice that ip± = 
on U\ A± to infer 

p 

Vs,p{w, Ax{U\A))< 2P~'^Vs,p{u, Ax{U\A)) + 2"-^ 



Eia(C/) 'Pi{x){u{x) - Zi) 



-dxdy 



Ax(c/\A) \x-y\^+-P 

We fix i e Z" and define A'^ = {Ai x (U \ 'A)) D Ap, then by using \Lpi{x) - (pi{y)\ < 2m'^\x - y\/p, 
with ipi{y) = for y £ U \ A, Fubini theorem and a direct integration yield 

Mii^..* < (2„.»)v- (3.17) 

/ /o 2nc -V f j f \u{x)~Zi\P nuJn{2m'^)P /" , . , 

<(2m^)Pp P I dx I , , I / \^ dy ^ r^p ^ / M(a;) - Zi rdx. 

Let now A'/ = (Ai x ([/ \ A)) \ Ap, then we argue as above using \(pi{x) — fi{y)\ < 1 to get again by 
a direct integration 

l^.(.)(u(.) -zor < ^ K-)--;^.,, 

< / / = ^P-' I H-) - -^A'dx. (3.18) 

By taking into account (3.16)-(3.18) we deduce 

Vs^p{w,Ax {U\A)) <2P-^Vs^p{u,Ax {U\A)) + c'm^Pp~'P ^ f \u{x) - Zi\Pdx. (3.19) 

By collecting (3.9), (3.11), (3.15), and (3.19) we infer w e W'-p{U) and estimate (3.8) for U xU. 
For any measurable subset E in U x U, being w = w on ?7 \ ^, we have 

\Vs,p{w, E) - Vs,p{u, E)\ = \Vs,p{w, En{UxA))- Vs,p{u, E n {U x A))\ 

^ I^IV-p(A) + T^sA^^ {U\A)xA)+ Vs,p{u, U X A). 
Eventually, (3.8) follows at once by taking into account (3.10)-(3.19). □ 
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By scaling Lemma 3.8 we establish a joining lemma for fractional type energies. 
Before starting the proof we fix some notation: having fixed m e N and set X, := T\.{U), for all 
i e Ij and /i e N let 

B}''* = {x e R" : \x- < Cf'' {x £ R" : m-^'^-V, < |x - xj| < m^^'^-Vj}. 

Clearly, we have C'/' C Bj^'' \ sj'''^^ C V/. 

Lemma 3.9. Let (uj)jgN &e converging to u in U'{U) with sup^ |Mj|i4/s,p([/) < +oo. With fixed m, 
€ N, for every j G N f/iere exists hj G {1, . . . , N} and a function Wj € W^'^lU) such that 

wj=u, onU\U,ei,(B^"' \B';''''^'), (3.20) 

w,{x)^{u,)^^..^ onCf'^', (3.21) 

for some c = c{n,p, s, m) > /loZds /or every measurable set E in U x U 

\Vs,piu,,E)-VsAw,,E)\ < (3.22) 

and the sequences (u'j)jgN; (CjOjgn, wii/i Cj •= 'l2i(EX -{u)i''^i) ^^•''3 Xv^ ! converge to u in LP{U). 
In addition, if Uj € L°°{U) 

WwjWl^^u) < \\uj\\l'^{u)- (3.23) 

Proof. Given m, N E N, then for every j £ N and He {1, . • . ,N} fixed, apply Lemma 3.8 with 
(A')i := C)-'', A'l := B'j''' = {uj)c^,h, i e X^. Take note that p = m-^'^r^-. if w]:^ 
denotes the resulting function and ~ Uigi^. then for some constant c — c{n,p, s) and for any 
measurable set E in U x U hy (3.8) it holds 

\Vs,p{uj,E)^Vs,p{wj,E)\<cVs.p{uj,U X A'') + cm'^P [ — ] V / \uj - {uj)(...HY'dx. 

This estimate, together with the scaled Poincare-Wirtinger inequality (2.8) with r = m~^^rj^ gives 

\Vs,p{uj,E)-Vs,p{wj,E)\ < c (Vs,p{u,,Ux A'') + |uj|^V=,p(a'-)) ^ cV,,p{u,,U x A''), (3.24) 

fro some c = c{n,p, s, m) > 0. By summing up and averaging on h, being the A'"s disjoint, we find 
hj G {1, . . . , A^} such that 

V,^p{uj, U X < j^V,^p{uj,U X UhA^). (3.25) 

Set Wj := w^j''^\ then (3.20) and (3.21) are satisfied by construction, and moreover (3.24) and (3.25) 
imply (3.22). 

To prove that (u)j)jgN converges to u in LP{U) wc use (2.8), with r = rji^^^rj, and the very 

definition of Wj as convex combination of Uj and the mean value {uj) i.hj on _B^'''^ \ B'j^'^^ to get 

j 
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where c = c{n,p, s, m) > 0. 

Eventually, let us show the convergence of (Cj)jeN to u in L'p{U). To this aim we prove that 
(Cj — Uj)jeN is infinitesimal in L^iU). Fix any number Af with sup^ Rj/rj < M < +00 (see (3.1)), 
we claim that for some constant c = c{n,p, s, m, M, N) > we have 

E ^ KOc'-"^ llLp(y.i) ^ crj^|uj|^,,p(y). (3.26) 

Given this for granted the conclusion is a straightforward consequence of the definition of Q, of (3.26), 
of (2.4)3 and of the equi-integrability of (|uj|'')jgN, i-e. 

lICi — '"■j\\^LP{U) ^ E ^ (^jOp'-'^i llLP(y^i) + I1^jIIlp(c/\ux. V/)- 

To prove (3.26) we use (2.8) and the fact that the balls i?A/rj (xj) have (uniformly) finite overlapping. 
More precisely, the inclusions V^- C Bii^{Xj) C Bmt ji^]) and (2.8) applied with r = rj give for some 
c = c{n,p, s, m, M, N) > 

Moreover, since Aj is a Delone set and by definition rj < r/y. , an elementary counting argument 
implies supigz„ #{k e Z" : Ba/,^ (xj) n Bmv, (x^) ^ 0} < (2M + 1)" for all 3 G N. 
Finally, (3.23) follows by construction. □ 

3.3. Proof of the F- convergence. We establish the F-convergence result in Theorem 3.3. It will 
be a consequence of Propositions 3.10, 3.11, below. We start with the lower bound inequality. 

Proposition 3.10. For every Uj u in U'{U) we have 

liminf Fj{uj) > J-{u). 
j 

Proof. Fix N E ISi , 6 > 0, and set m ~ li-/S\ G N, [-J denoting the integer part function. Consider 
the sequence (u'j)jgN provided by Lemma 3.9. We do not highlight its dependence on 6, N for the 
sake of notational convenience. We remark that whatever the choice of S and N is, we have that 
(wj)jizT^ converges to u in U^{U) and that for some c ~ c(n,p, s, 6) it holds 

(1 + —) liminf Fj(uj) > liminf Fj(wj). (3.27) 

\ N / j j 

Furthermore, we note that for j sufficiently big Uigi^ {V^^ x Vj-) C As, and thus 

liminf F,(7i;,) > liminf | / ^M^ifiM^d^dy + y \yj \P 
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>/ l^^(a^) "yi^ rfa^d^ + liniinf y M^^,,,,,, (3.28) 



thanks to Fatou lemma. We claim that 

Immf > e (cap,,,(r) - e,) ^ Hx)^ P{x)dx, (3.29) 

with £5 > infinitesimal as i5 ^ 0^. Given this for granted, by (3.27) inequality (3.28) rewrites as 

(l + |^)limmfi^,(u,)> 1^^^^^ '"11" |n£r + 9 (c^p^jT) - es) j^\u{x)\^ I3{x)dx. (3.30) 

The thesis then follows by passing to the limit first as iV — > +oo and then as (5 ^ 0+ in (3.30). 

To conclude we are left with proving (3.29). We keep the notation of Lemma 3.9, and further set 
B] := {x e R" : \x - xj| < m^^^/^.+i)^^}^ fo^, i e Xj. Take note that B] C Vf. We have 



',C(v;.) > mf : w € W^^'^'(R"), w = iu,)^., on C^'^ w = q.e. on Tj j (3.31) 



inf <{ : w e W'^'f(R"), w = Oon C;■'^ u; = (t.,)^,., q.e. on 



= reap,, (l-, . Ar-|(..),;.S reap,., (^T, 5-^; ^^^3^%) ■ 

The last equality is justified by an elementary translation and scaling argument. Thanks to (2.17) in 
Lemma 2.12, by recalling that hj G {1, . . . , N}, we get the following estimate 

cap, T,B ; ^ra^ > cap,.p(r) - ^-p— C^.p Bi,B . 



Hence, if A e A{U) is such that A CC ?7, for j sufficiently big we infer 

E k.r„..,.(^.) > (^cap,.,(r)-^^^^^C.,p I^MxV^,ix)dx, 



where ^j(x) := X^iei -^j *^('C"(Vj^)) ^XyK^^) ^.nd Q is defined in Lemma 3.9. Take note that by 
(2.4)2 we have 



^,{x)dx~yr''mA'}) 



A' 



for any A' G A{U). Testing the inequality above for all cubes in U with sides parallel to the coordinate 
axes, centers and vertices with rational coordinates yields j 9 (3 weak* L°°{U) by (3.2) and (3.3). 

By taking this into account, the thesis follows at once by the convergence of relative capacities 
to the global one proved in (2.16) of Lemma 2.12, the strong convergence of (Cj)jeN to u in LP{U) 
established in Lemma 3.9, and eventually by letting A increase to U. □ 
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In the next proposition wc prove that the lower bound established in Proposition 3.10 is attained. In 
doing that the main diffieulty is to show that the inequalities in (3.28) are optimal. Due to the insight 
provided by Proposition 3.10 we show that the capacitary contribution is concentrated along the 
diagonal set A and is due to short range interactions. Instead, long range interactions are responsible 
for the non-local term in the limit. 

Proposition 3.11. For every u G U'{U) there exists a sequence {uj)j^TSi such that uj u in U'[U) 
and 

lim sup (wj ) < J-{u). 

j 

Proof. We may assume u € W'^'°°{U) by a standard density argument and the lower semicontinuity 
of P- lim sup Tj . 

With fixed € N, let ^jv G W^'*-p(R") be such that = on R" \ B^, > 1 cap, ^ q.e. on T 
and ICwl^s.p(R^n-) < Cs,p{T, Bn) + 1/N, and let S C^{Bn) be any function such that = 1 on Br^, 
(LipC)P < 2 and < C < 1- 

Let (wj)jgN be the sequence obtained from u by applying Lemma 3.9 with m — 2. We keep the 
notation introduced there and further set 

:= {a; G R" : |x - xj| < 2^^'^'rj}, := u^i,h^ for every i G X,-, 



:= Bnx^ (xj) n U for every i e ^j, := ^a, (C/), 
C/,:=L/\((Ux,Sl)u(u^Bj)^ 



Then, recalling that Xj — rj^'" define 



Uj{x) := < 



Wj{x) Uj 

^l-C^(^))"} B^^tel, (3.32) 

For the sake of notational simplicity we have not highlighted the dependence of the sequence (itj)jgN 
on the parameter G N. Clearly, (Mj)jgN converges strongly to u in LP{U), and moreover it satisfies 
the obstacle condition by construction. The rest of the proof is devoted to show that Uj G W'PlU) 
with 

limsupJFj(itj) < T{u) + es + eat. 



3 



where — > 0+ as (5 ^ 0+ and ejv ^ 0+ as ^ +oo. 

A first reduction can be done by computing the energy of Uj only on a neighborhood of the diagonal 
A. Indeed, Lebesgue dominated convergence and the stated convergence of (uj)jgN to u in LP{U) 
imply 

liml?,,p(uj, iUxU)\ As) = V.^piu, iUxU)\ As). 

3 
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In addition, since Uj = Wj on Uj by (3.22) in Lemma 3.9 we have for some constant c = c{n,p, s) 

\imsupVs,p{uj, {UjXUj)nAs) < limsupl?,,p(wj, {UxU)nAs) < (^1 + ^) iUxU)nAs) = eg. 

(3.33) 

The conclusion then follows provided we show that 

limsup (2?.,p(wj, iUx{U\ TJj)) n As) + ^^^^(wj, {{U \ Uj) x Uj) n As)) 
j 

<ecap,p{T) [ \u{x)\Pf]{x)dx + eN + es. (3.34) 



In order to prove this we introduce the following splitting of the left hand side above: 

I?.,p(u„ {Ux{U\ U,)) n A,) < ^ + E T^sM,,B^ X B^) 



iGl,- 



{(i,k)el|:0<|xj-x5|<5} 



+2 E X [/,) n Ai) + 2 E ^^^P (^J' (^1 ^ ^j) =• ^1 + • • • + ^: 

1G.#3 (i.k)elj x,J?j 



Next we estimate separately each term /j*, h E {1, . . . ,6}. Since the computations below are quite 
involved, we divide our argument into several steps to provide a proof as clear as possible. Take 
note that all the constants c appearing in the rest of the proof depend only on n, p, s, hence this 
dependence will no longer be indicated. 



Step 1. Estimate of I. 



1 . 



limsup/j < 6'(cap,p(T) + ejv) /" \u{x)\p I3{x)dx. (3.35) 
' Ju 



A change of variables yields 



(A7^(Bi-xj)) 



\Qix)\P^,ix)dx, 

u 



where j{x) = X^iei ■^j {i^^{Vj)) '^xyi{x) and Q is defined in Lemma 3.9. Arguing as in Propo- 
sition 3.10 and by Lemma 2.12 we conclude (3.35). 



Step 2. Estimate of I'j . 



limsup/j^ < es- (3.36) 



Take note that by the very definition of uj in (3.32) for any {x,y) G i? j x Bj, i k and i,k e X, , 
we get 

M^) - uM - [u] - u)) - {Xj\x - xj)) u) + {Xj\y - xj)) u^;. 
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Hence, we can bound ij as follows 



ielj {kGlj:0<|xi-x5|<5} 



Since |xj — ^jl/^ !i k — y| < 2|x^ — x^| for any {x,y) £ Bj x B^, i,k e X,- with i 7^ k, we infer 
J — u^l < 2Lip(w)|xJ — x^l being u £ Lip(J7), and so we deduce 



1"}-"^!" ........ 



2n 



< cUpP(u)— rf-rr— T^- (3-37) 

|a; - ^- ^ ^ |x^ - x^|"+(^-i)p ^ ^ 

J j J J 

To go on further wc notice that for every fixed i € Xj we have 

{k e X, : < Ixj - x^U <S}(Z uitf^^J {k € I, : hr, < jx) ~x^\^<ih + l)r,}, 

where [tj denotes the integer part of t. The latter inclusion together with (2.4)i, (2.5) and (3.37) 
entail 

[S/r.,] n-(s-l)P 

^ ^Lip^(-) EE E 

ielj h=2 {keIj:/irj<|xi-x5U<(/i+l)rj} 

(2.4)„(2.5) L^A^J . , 

< cLip^'(^) ^ < cLip''(u)<5^i-^)^ (3.38) 

h=2 

In the last inequality we used that X]/!=2 < (-^^ ''')/(~7); for ^-ny 7 < and M G N. 

To deal with I^'^ we use a similar argument. Indeed, for every i S X, we have 



E 



V /" 1 ^ V- (2.5) 1 



Thus, being ^Ar(A ■ ^(- — x^)) supported in i?^, a change of variables yields 



3' 

. Ti-sj; 



ir<-\\<^^nn^''''>^7rT--^'U^\\^^^^^ (3-39) 



Clearly, (3.38) and (3.39) imply (3.36). 
Step 3. Estimate of : 



limsup/^ < e^- + ejy. (3.40) 



Being Uj = Wj on Uj and spt(^Ar(Aj ^(- — xj))) C Bj, we find 

l! < chll^^(^) E T^s,p(iN{Xj\--^)lB^ X ([/\ B))) 
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7"3,1 I 7"3,2 I 7-3,3 

Note that by a change of variables the integral rewrites as 

< c\\u\\l^^^^X]-'^r--V.,,p (^Cn,B^ X I^R" \ = CAT, (3.41) 

by Remark 2.14. To deal with the term Ij'^ we first integrate out y and then use Hardy inequality: 

- ^ ^ L disl-^(x,9i^^) ^" - ^ ^ l-^ l'v.=..(i.)) ^ X ^) n = (3-42) 

Finally, for what /^'"^ is concerned we have 

< cP,,p(u.j, (t/ X C/) n As) 65. (3.43) 
By collecting (3.41)-(3.43) we infer (3.40). 
Step 4- Estimate of ij: 

limsup/j' < £5. (3.44) 

j 

The very definition of Uj in (3.32) gives for any {x,y) E Bj x Bp i,k G J^j 

u,{x) - u,{y) = (1 - C {Xf{x - xj))) w,{x) - (1 - C {\-j\y - x^))) w,{y) 

= (1 - C (A-^(x ~ xj))) (z«,(x) - w,{y)) + (C (A-^(x - x^^)) - C {)^]\y x^))) ^M- 
Distinguishing the couples of the form (i, i), i e from the others we bound as follows 



p 



A change of variables and (2.4)2 yield 

(.3.33) 

Ip'+Ip' <cV,A^,,iUxU)nAs)+cC^^{{dU)R^)\Cf^^,,^^^^ < es + cr^{idU)R^). (3.45) 
To deal with the term ij'^ first note that 

then integrate out y, scale back the x variable, and finally use Hardy inequality taking into account 
that c e C^{Bn). 

<c||z.||,^(^,^|^^y^^^ dist-(x,9i?j) ^|^^y^^dist-(.,9B^)^" 
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by (2.4)2. In conclusion, (3.45) and (3.46) give (3.44). 
Step 5. Estimate of Ij : 

limsup/| < £5. (3.47) 

j 

The computations are similar to the previous step once one notices that for any {x,y) (z Bj x Uj, 
i G it holds 

Uj{x) - Uj{y) = (1 - C {>'J^ix - x}))) Wj{x) - Wj{y) = (wjix) - Wj{y)) - ( {>^J^ix - x))) Wj{x). 

Hence, we bound /| by the sum of two terms, the first analogous to Ij'^ and the second to ij'^. Thus, 
(3.47) follows. 

Step 6. Estimate of : 



J 

1 mi ciTr\ I . 

3 



limsup/f<e5. (3.48) 



For {x,y) £ Bj X Bj, i e Ij and k G J^j, we write 

Uj{x) - Uj{y) = (1 - Cat {\~^{x - x^"))) [u] - Wj{x)) 

+ (1 - (A7^(x - 4))) w,{x) - (1 - C {Xf{y - 4))) w,{y). 

Thus, we infer 

^^^^ 2^ /. \r y\n+sp dxdy + c\\ur ^ / ' ' dxdy 

^ II IIP \^ f \^n{X-\x~^)))\p 

Clearly, /|'^ can be estimated as Ij'^, Ij'^ as ij'^, and /|'^ as In conclusion, (3.48) follows. 
Step 7: Conclusion. By collecting Step 1 - Step 6 we infer 

limsupJFj(itj) < .F(u) + es + eN, 

j 

with the two terms on the rhs above infinitesimal as i5 — s- 0^ and as — s- +oo, respectively. □ 

3.4. Generalizations. Anisotropic and homogeneous variations of the fractional semi-norm can be 
treated essentially in the same way. Consider a ^"-measurable kernel K : R" \ {0} (0, +oo) such 
that for all z G R" \ {0} and for some constant a > 1 it holds 

A'(tz) = t-("+"f)X(z) t>0, a~^\z\-^''+''P'> < K{z) <a\z\-^''+'Pl (3.49) 

Define JC : ^^^^^(R") x y^(R") ^ [0, +oo) as 

JC[u,A):= I K{x-y)\u[x) -u{y)\P dxdy, (3.50) 

J AxA 
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dropping the dependence on the set of integration in case A = R". Ah the relevant quantities 
introduced in the preceding sections have analogous counterparts simply replacing the kernel | • 
with K. For instance, the locality defect associated to the energy K. is given by 

Vk.{u,E) := I K{x~y)\u{x)~u{y)\Pdxdy, 

J E 

for any E C- U x U "-measurable. We point out that in general, V/ciu, A x B) ^ Vjciu, B x A), 
A and B being measurable subsets of R". Nevertheless, a splitting formula similar to (2.15) holds. 

The only relevant changes are in the proof of Lemma 2.12 in which we exploited the invariance of 
the kernel | • under rotations to establish (2.16). As already noticed before the statement of 

Lemma 2.12. and as it turns out from the proof of Propositions 3.10, 3.11, in the current deterministic 
setting it is sufficient to prove pointwise convergence of the relative capacities to infer the lower 
bound estimate. This is the content of the next lemma. The argument below does not give uniform 
convergence. 

Lemma 3.12. Let r > and define for every T C R" 

CK.iT, Br) :=inf{/C(w) : we VF^'P(R"), w = on R" \ B^, w > 1 q.e. on T} . 

Then 

lim CKiT,Br)^cap,^iT). (3.51) 
Moreover, there exists a constant c — c(n^ s^p^ a) such that for all < r < R 

cap^(r) -cap^(r,i?H;r) < ^-^I1—Ck{T, B,.), (3.52) 

where capj^(r, B^; := inf {/C(u;, 5^) : w £ W'P(R"), w = on W \ Br, w > 1 q.e. on T}. 

Proof. Estimate (3.52) can be derived as we did for (2.17) thanks to (3.49)2. 
It is clear that cap^iT) < Cic{T, Br) < Ck.{T, Br) for all < i? < r, so that 

lim CK{T,Br) > cap^(r). 

r — !-4-c>o 

Fix r > such that T C Br/2 and let u e W'^p{W) be such that < m < 1 £" a.e. on R", u > 1 q.e. 
on T . Given (pr a cut-off function between Br/2 and B^r/A set Wr :~ iprU. By construction Wr = u on 
Br/2, = on R" \ B:ir/i and Wr > 1 q.e. on T. We claim that Wr G M^*'P(R") with 

IC{Wr,Br) < K(u,Br/2) (3.53) 



+c(n,p, s, a) Vk{u, Br/2 x (B, \ Br/2)) + Vk{u, (Br \ Br/2) x Br/2) + r''"' / _ \u\Pdx . 

\ JB^\Br/2 J 

To this aim we bound the energy of Wr on Br by 

IC{Wr, Br) < IC{U, Br/2) +VK{Wr, Br/2 X (B, \ B./a)) + Pyc (u^r, (S, \ B./a) X Br/2). (3.54) 
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We estimate the first of the two locality defect terms in (3.54), the argument for the second being 
analogous. In doing that we can follow the argument used in Lemma 3.8 for the locality defect term 
(see (3.17)-(3.19)). Then, by (3.49)2 and since < ipr < 1, we infer 

< 2P Vk{u, Br/2 X {Br \ Br/2) + « / _ 1 d^^V 



< c{n,p, s, a) Vk{u, Br/2 x {Br \ Br/2) + r-^" / _ \u\Pdy . 

\ Jb,\b,^2 J 

Formula (3.53) then follows at once. 

Furthermore, since Wr — on -83^/4, (3.49)2 and (A. 2) in Lemma A.l yield 

V!c{wr,Br X (R" \ B^)) < c{n, p, s)a r-P' f \u\Pdx, 

J Br 

and thus we conclude Wr E VF^'P(R"). 

Eventually, since lim^ X>k;(w, R" x (R"\Br/2)) = lim^ 'Dic{u, {W\Br/2) x R") = and u G Lp(R"), 
from (3.53) we conclude 

liinCic{T, Br) < limsup/C('u;r, -Br) l£ K,{u). 

r r 

Taking the infimum on all admissible functions u we conclude. □ 

With fixed U G ^(R"), consider ICj : LP{U) [0, +00] given by 

I JC{u, U) if u e Ty"^P(C/), u = cap, q.e. on Tj n U 
lCj{u)=l ' (3.55) 

+00 otherwise. 

Theorem 3.13. Let U G ^(R") be bounded and connected with Lipschitz regular boundary. 

Then the sequence (/Cj)jgN T-converges in the LP{U) topology to J(f : LP{U) [0, +00] defined by 



jr(u) = /C(u) + 6icap;c(7') / \u{x)\P(3{x)dx (3.56) 

Ju 

if u G W^'P{U), +00 otherwise in U'{U), where 

cap^(r) inf {/C(w) : w G W^^'P(R"), w > 1 q.e. on T} . 

Proof. The proof is the same of Theorem 3.3 a part from the necessary changes to the various technical 
lemmas preceding Theorem 3.3. We indicate how these preliminaries must be appropriately restated. 

We have shown the (pointwise) convergence of relative capacities in Lemma 3.12. Changing the 
relevant quantities according to the substitution of the kernel | • with K, Lemmas 3.8 and 3.9 

have analogous statements since the splitting formula (2.15) does. 

In the proof of Proposition 3.10 we use the homogeneity of the kernel K (see (3.49)i) for the scaling 
argument in (3.31) leading to the analogue of formula (3.29). This is the reason why we ask for (3.49)i. 
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Proposition 3.11 needs few changes: in the definition of uj choose ^jv to be a (l/A'^)-niinimizer 
of Cfc{T, Bn), the estimate of ij follows straightforward. For what the terms /j', h £ {2, ... ,6} are 
concerned it suffices to take note that by condition (3.49)2 the locality defect satisfies a^^Vg.piu, E) < 
VjciujE) < a'Ds,p{u,E). Hence, we can follow exactly Steps 2-6 to conclude. □ 

Remark 3.14. Dropping the homogeneity assumption (3.49)i on the kernel K the T-limit might not 
exist. This fact had already been noticed in the local case (see [')]). Though, in such a framework 
abstract integral representation and compactness arguments for local functional imply T-convergence 
up to subsequences. 

4. Random Settings 

In this section we extend our analysis to two different random settings. First, we deal with obstacles 
having random sizes and shapes located on points of a periodic lattice as introduced by Caffarelli 
and Mellet [18], [19] (see also [32]). In particular, we provide a self-contained proof of the results in 
[19], [32] avoiding extension techniques. Second, we consider random homothetics copies of a given 
obstacle set placed on random Delone sets of points following the approach by Blanc, Le Bris and 
Lions to define the energy of microscopic stochastic lattices [13] and to study some variants of the 
usual stochastic homogenization theory [14]. 

We have not been able to work out a unified approach for the two frameworks described above. 
The main issue for this being related to the interplay between the weighted version of the pointwise 
ergodic theorem in Theorem 4.2 below and stationarity for random Delone sets of points (see (4.16)). 

In both cases we are given a probability space (ri,^,P) such that the group Z" acts on fl via 
measure-preserving transformations ti : — > 51. The cr-subalgebra of ^ of the invariant sets of the 
Ti's, i.e. O C n such that tiO = O for all i e Z", is denoted by J^. Recall that (Ti)igZ" is said to be 
ergodic if y is trivial, i.e. O £ J satisfies cither P(0) = or P(0) = 1. 

In the sequel we keep the notation introduced in Section 3 highlighting the dependence of relevant 
quantities on lo when needed. 

4.1. Obstacles with Random sizes and shapes. In this subsection we deal with the case of 
obstacles with random sizes and shapes located on points of a lattice. We restrict to the standard cubic 
one only for the sake of simplicity (see Remark 4.3 below for extensions). More precisely, let = £jZ" 
with (£j)jgN a positive infinitesimal sequence, then we have x^ = £ji, = £j{l + [—1/2,1/2]"), 
rA^ = ej/2 and Ra^ = In addition, if = ej/2 then 6 = /:"(;7)/2" and /? ee l/C"{U). Thus, 

to simplify the presentation we will use the scaling parameter ej instead of rj , and the more intuitive 
notation for Vj-. 

Let us now fix the assumptions on the distribution of obstacles originally introduced by Caffarelli 
and Mellet [IS], [19]. For all u; e fl and j G N the obstacle set Tj{uj) C R" is given by T,(tj) := 
^i£Z"Tj{uj), where the sets Tj{uj) C Q^, satisfy the following conditions: 
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(01) . Capacitary Scaling: There exists a process 7 : Z" x ^ [0, +00) such that for aU i G Z" and 

cap,^p(r;H)=£;'7(i,w)- 
Moreover, for some 70 > we have for aU i e Z" and P a.s. ui G ft 

7(i,w)<7o. (4.1) 

(02) . Stationarity of the Process: The process 7 : Z" x ^ [0, +00) is stationary w.r.to the family 

(Ti)igz", i-e. for all i,k e Z" and uj E ft 

7(i + k,cj) =7(i,Tkw). (4.2) 

(03) . Strong Separation: There exists a positive infinitesimal sequence ((5j)jgN, with 6j = o{ej) and 

e]^^ = 0{Sj), such that Tf{uj) C xj + dj{Q) - xj) for all i e Z", w e rj. 

Take note that by (02) we have 7(1, lu) ~ 7(0, t-±lu), hence the random variables 7(1, uj) are identically 
distributed. The common value of their expectations is denoted by £[7], i.e. 

E[7] -.^ / 7(0,w)dP. 
Jn 

Moreover. E[7, J^] denotes the conditional expectation of the process, i.e. the unique .^-measurable 
function in L^{il,P) such that for every set O G 



7(0,Lj)(iP= / E['y,J^]{Lu)d¥. 
' Jo 

In the sequel we analyze the asymptotics of the energies ICj : LP{U) x ^ [0, +00] given by 

\JC{u) if u e W'^P{U), u^O cap q.e. on Tj{oj) n U 
JCj{u,uj)^< ^ ' (4.3) 

+00 otherwise, 
where /C is the functional in (3.50). 

Theorem 4.1. Let U e ^(R") be bounded and connected with Lipschitz regular boundary, and assume 
that the kernel K satisfies (3.49) and is invariant under rotations. 

If (Ol)-(OS) hold true, then there exists a set 17' C of full probability such that for all uj G 17' the 
sequence (/Cj(-, w)) T-converges in the LP{U) topology to J(f : LP{U) x fl^ [0, +00] defined by 

jr(u,w) =/C(m)+E[7,^] /" \u{x)\Pdx (4.4) 



ifu e W'^P{U), +00 otherwise in LP{U). 
If in addition (Ti)igZ" is ergodic then the T-limit is deterministic, i.e. £[7,^] ~ £[7] P a.s. in Q. 

The proof of Theorem 4.1 builds upon Theorem 3.13 and upon a weighted ergodic theorem estab- 
lished in [32, Theorem 4.1]. We give the proof of the latter result for the sake of convenience. 
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Theorem 4.2. Let 7 be satisfying (02), then for every hounded open set V C R" with C"'{dV) = 
we have P a.s. in 

lim— ^ J2 7(i,c^)=E[7,^], (4.5) 



^jix^u):^ J2 7(i,'^)XQj(^)^E[7,^] weak*L-(F). (4.6) 

±eij{V) 

Proof Define the operators Ti : L^{n,P) L\n,F) by Ti(/) -.^ f o n for every i e Z". The group 
property of (Ti)igZ" imphes that = {Ti}^^^^ is a multiparameter semigroup generated by the 
commuting isometries T^^ for r e {1, . . . , 71}, being {ei, . . . , e„} the canonical basis of R". 

We define a process F on bounded Borel sets V of R" with values in L°°{n,¥). Set Qi := 
[-1/2, 1/2]" and let F be as follows 

Fv{uj) := ^ 7(1, w), 
{ieZ":i+QiCV} 

with the convention that Fv{lo) := if the set of summation is empty. It is clear that F is additive, 
that is it satisfies 

(i) F is stationary: Ti o Fy ^ Fy+i for ah i e Z"; 

(ii) FV1UV2 = Fv, + Fy^, for disjoint Vi, V2; 

(iii) the random variables Fy are integrable; and 

(iv) the spatial constant of the process 7(F) := inf{j~" FjQ-^dF} G (0, 00). 

Indeed, (i) and (iv) follow by (02) (actually 7(F) — IE[7]), (ii) by the very definition of F, and (iii) 
by the positivity of 7 and (4.1). It then follows from [■>(), Theorem 2] and [3(), Remark (b) p. 294] that 
there exists / e L^{n,F) such that for all V G 6(R") bounded with C"{dV) and P a.s. in n 

lim j-^F,v{uj)^C^{V)7{uj). (4.7) 

Stationarity and boundedness of 7 together with (4.7) yield that for any O e =^ we have 

7{uj)d¥^ \im / FjQ,{oj)dF= f 7(0,c^)dP, 



so that the limit / is actually given by E[7, J^] (see for instance [■!'), Theorem 2.3 page 203]). 

In particular, (4.7) still holds by substituting (j)jeN with any diverging sequence (aj)jgN Q N- Take 
Oj = [l/^iJ (W stands for the integer part of t), then notice that for every S > and j sufficiently 
big we have 



7(i, w) - i^a,y(w) 

iGl,(V) 



< 7oa-"# {ljiV)A{± e Z" : aj\± + Q,) C V}) < ^oC\{dV)s), 



here A denotes the symmetric difference between the relevant sets. Since e"#X,(l/) £"(y) and 
QjEj — > 1 we infer (4.5). 
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Eventually, in order to prove (4.6) consider the family £S. of all open cubes in R" with sides parallel 
to the coordinate axes, and with center and vertices having rational coordinates. To show the claimed 
weak* convergence it suffices to check that P a.s. in it holds 

lim / ^,{x,w)xQ{x)dC' ^C\Q)Wa.A 
3 Jn 

for any Q E ^ with Q C X^, 1/ as in the statement above. We have 



Jo 



< 



27o/:"(Q\Uiei^.(Q)Qj) 



and thus (4.5) and the denumerability of ^ yield that the rhs above is infinitesimal P a.s. in il. 



□ 



Proof (of Theorem 4-6)- We define fi' as any subset of full probability in for which Theorem 4.2 
holds true for U . Then, we follow the lines of the proof of Theorem 3.f 3 pointing out only the necessary 
changes. 

First, we note that by assumption (03), in particular condition Sj = o{ej), we can still apply 
Lemma 3.9 in this framework. Hence, to get the lower bound inequality we argue as in Proposition 3.10, 
we need only to substitute (3.29) suitably. Formula (3.31) is replaced by 



w,\L..,,y., > e"|(«,) ...reap,,, XjHniu;)-x)),B_ 



>(7(i,c.)-e„) s]\{u,) \P, 

3 

where > is infinitesimal as to ^ +oo. To infer the last inequality we have taken into account 
(01) and the uniform convergence of the relative capacities established in (2.17) of Lemma 2.12. This 
is guaranteed by condition ^^^"--■'p — 0{5j) in (03), which ensures that the rescaled obstacle sets 
XJ^{Tj{lj) — x^) are equi-bounded. 

By summing up all the contributions, by Theorem 4.2 and by recalling that the sequence (Cj)jgN 
defined in Lemma 3.9 converges strongly to u in U'{U) we infer 

liminf y > liminf / ,{x,uj)\Cj\Pdx > (E[7, J^] - e™) / \u{x)\Pdx, 

for all A e A{U) with A CC U . By increasing AtoU and by letting to +oo we conclude. 

The upper bound inequality is established as in Proposition 3.11 by substituting S,n in the definition 
of Uj with Cj'^ y ^ (l/Af)-minimizcr of Ck;(AJ^ (Tj^((jj) — Xj),Bn). The uniform convergence of those 
relative capacities is guaranteed by the analogue of (2.16) in Lemma 2.12. This is the reason why we 
suppose K to be invariant under rotations. 
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The estimate of /j then foUows straightforward. For what the terms /j*. h G {2, . . . , 6} are conecrned, 
take note that sinee a^^Vg.piu, E) < Vjciii, E) < a'Ds^p{u, E) (see (8.49)2), we can follow exactly Steps 
2-6 to conclude. □ 

Remark 4.3. If A is a generic periodic n-dimensional lattice in R" we can argue analogously. By 
definition the points of A belong to the orbit of a Z-module generated by n linearly independent vectors, 
and the Voronoi' cells turn out to be congruent polytopes (see [42, Chapter 2]). Clearly, the obstacle set 
Tj[ijj) and stationarity for the process 'y{-,iu) have to be defined according to the group of translations 
associated to vectors in the mentioned Ti-module which leave A invariant. 

4.2. Random Delone set of points. According to Blanc, Le Bris and Lions (see [13], [14]) a random 
Delone set is a random variable A: Vl ^ (R")'^ satisfying Definition 2.2 P a.s. in Vl. Then, both ta 
and i?A : r2 ^ R turn out to be random variables: ta because of its very definition (2.2)i, and i?A 
since it can be characterized as Ra{') ~ supq„ dist(.T, A(-)). 

We will deal with sequences Aj : Q (R")^ of random Delone sets fulfilling conditions analogous 
to (3.1)-(3.3) (see below for relevant examples), that is for P a.e. lj e it holds 

lim||rj||L^(o.p) = 0, (1 <) linisup ||i?j7rj||i^.(o.p) < +00, (4.8) 
^ j 

lim#(A,(cj) n C/)r"(tj) = 9{uj) e (0,+oo), (4.9) 
j 

:= M^.(\r.u) E <5.jM(-)^M-,^) :-/3(-,^K"LC/ w*-Ct{U), (4.10) 

for some B{U) ® !P measurable function /3 such that /3(-,ct;) G L^{U, [0, +00]) and ||/3(-, cj)||ii(^r') = 1. 

Remark 4.4. In (4.10) the set of indices Aj{uj) D U sushsitutes Xj{U, lu) to ensure the measurability 
of 6, (3. For, the measurability of ff:{A{Lu) D U) follows easily from that of the random Delone set A. 
Despite this, since < #{Aj{lo) HU) - #Jj{U,lu) < if.yj{U,uj), by (4.8) we have for P a.e. uj e n 

.^.i^ . E w*-Ci,iU), l\inifI,iU,iu)r^ico) = 0{io). (4.11) 

Remark 4.5. Contrary to the deterministic setting we do not know whether conditions (4.9), (4.10) 
are satisfied up to the subsequences or not. Nevertheless, in Examples 4-9, 4. 10 we show some sets of 
points satisfying all the conditions listed above. The only piece of information we extract from (4.8) 
is that for a subsequence ( not relabeled in what follows ) we have the convergence 

{h{-,^),^)c,{U),v{u)~^ {l^{-,^),v)c,{u),v(U) w*-L°^{n,¥) for all^eCbiU). (4.12) 
Thus, condition (4.10) is equivalent to the strong convergence in L^{fl,P) for all ip e Cb{U) of the 
(sub) sequences {{^.]{-,i^),^)ct(U),r{,U))- 

To establish (4.12) take note that fij is a Young Measure (see [S]j. More precisely, fij : {Q, ^,¥) 
'P{U), being V(U) the space of probability measures on U, is a measurable map, i.e. fj,j{A,u!) is 



APERIODIC FRACTIONAL OBSTACLE PROBLEMS 31 

measurable for all A £ B{U). By taking into account the uniform conditon (4.8) and by arguing as 
in Remark 3.2 we infer that the family (/ij)jgN is parametrized tight (see [8, Definition 3.3]J; given 
(5 > for a compact set Cs d U it holds 

{F(»fij){n xCs)^ [ Hj{Cs,uj)dF{uj) >1-S. 

Thus, by parametrized Prohorov theorem (see [n, Theorem 4.8] j (4.12) holds true up to a subsequence. 

Finally, in view of Remark 3.2 and (4.8) we note that it is not restrictive to suppose that the limit 
measure in (4.10) is absolutely continuous w.r.to C"\—U P a.s. infl. 

We fix a bounded set T C R" and define the obstacle set Tj{u!) := U±^z"Tj{u!) as the union of 
random rescaled and translated copies of T according to (3.5), i.e. 

T/(cj) := xiiuj) + \j{uj)T, where X,{uj) := (cc;)"/("-^p) . (4.13) 

The asymptotics of the energies ICj : LP{U) x Q, ^ [{), +oo] given by 

I K,{u) if u e W'PiU), u^Q cap, ,„ q.e. on TAuj) n U 
lCj{u,uj)^l (4.14) 

+00 otherwise. 

where /C is the functional in (3.50), is a straightforward generalization of Theorem 3.3. 

Theorem 4.6. Let U e ,4(R") be bounded and connected with Lipschitz regular boundary. 

Assume Aj : (R")^ is a sequence of random Delone sets satisfying (4.8)-(4.10) P a.s. in f2. 
Then, P a.s. in the sequence [K.j{-,Lo)) T-converges in the LP(U) topology to : LP{U) x fl ^ 
[0, +oo] defined by 

J(r{u,Lu)=K.{u)+e{Lu)capi^{T) I \u{x)\p p{x,uj) dx (4.15) 

Ju 

ifu £ W''P{U), +00 otherwise in LP{U). 

To ensure that the F-limit is deterministic we introduce a sequential version of stationarity for ran- 
dom lattices as defined by Blanc, Le Bris and Lions [13], [14]: there exists a positive and infinitesimal 
sequence ((5j)jgN such that for all i e Z" and P a.s. in il. 

Aj(TiLj) === Aj(tj) - (4.16) 

Corollary 4.7. // the assumptions of Theorem 4-6 hold true, if (Aj)jgN satisfies (4.16), and if the 
family (Ti)igZ" is ergodic, then ta^ is constant P a.s. in ft for all j G N, and there exists (3 G L^{U) 
such that for ¥ a.e. uj £ we have (3{-,uj) ~ f3 £" a.e. in U. 

In addition, if rj is constant P a.s. in f2 for all j G N, then 6 is constant P a.s. in 51. In particular, 
the T-limit in (4.15) is deterministic, i.e. P a.s. in Q, the functional in (4.15) takes the form 

Jf{u,u)=K.{u)+ec&Y>K,{T) [ \u{x)\Pf]{x)dx 

Ju 

for every u G W^'^iU), +oo otherwise in LP(U). 
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Proof. Wc show that r\. is invariant under the action of (Ti)igZ", this suffices to conclude since 
ergodicity imphes that the only invariant random variables are P a.s. equal to constants. We work 
with I j {A, Lu) in place of Aj{Lu)nA, since < #{Aj{uj)nA)-#Ij{A,uj) < #y^{A,Lu) for all A e A{U) 
as already pointed out for A = [/ in Remark 4.4. 

Given any i e Z", (4.16) yields Vj^{TiUj) = V^^{u)) — iSj for some ki S Z", which implies in turn 
that r\. (TiUj) — (lo), and thus r^j is equal to a constant P a.s. in il for every j G N. 

In addition, we have also that t±u!) = ^Tj{A+i5j,uj) for every A G A{U), then if = 

for any (5 > and j sufficiently big it holds ^Ij{A_s,bj) < ^Tj{A,T±u!) < ^Tj{As,uj) (see (2.1) for 
the definition of Ag and A^s)- In particular, wc infer 

m,{A,nuj) - #I,{A,uj)\ < #(Xj U ^j){As\Ato) V #(Jj U J^j){A\—s,u). (4.17) 

If e 17 is such that (4.10) holds, for any xq e U and r e (0, dist(a;o, 9?/)) we have by (2.4), (4.11)i 
and (4.17) 

In turn, from the latter equality we infer that if 

;5(a;o, w) := hmsup -r I3{x,uj)dx, 

then j3{xQ,bj) = f3{xo,T±uj) for all i G Z", G ?7 and for P a.c. cj G 51. Thus, ${xo,uj) is P a.s. equal 
to a constant for every xq € U; Lebesgue-Besicovitch differentiation theorem yields the conclusion. 

Eventually, suppose that rj is constant, then (4.17), C^{dU) = and (4.8) imply that 9{t±uj) ~ 6{lo) 
for all i G Z". In conclusion, 6 is equal to a constant P a.s. in Vl. □ 

We discuss some examples related to sets of points introduced in [I.?] and [14]. As before, (£j)jgN 
denotes a positive infinitesimal sequence. 

Example 4.8. Let us consider ensembles of points which are stationary perturbations of a standard 
periodic lattice. More precisely, given a random variable X : {Q,^,¥) —>■ R" define the family 
(Xi)igz" by Xi{Lu) := X{Ti{Lu)). By construction (Xi)±^z" is stationary w.r. to (Ti)i£Z"- Let A{lu) 
{x^(a')}igZ" where x.^ (uj) := ±+Xi(Lu), and Aj{uj) := ejA(w). A simple computation shows that (4.16) 
is satisfied with 5j = ej. Moreover, if M := sup^^ jl^i — A'H^ooj-fj p-j < 1, then A is a random Delone set 
with 1 — M < rA < Ra < I + M V a.s. in Cl. In conclusion, (Aj)jgN satisfies both (4.8) and (4.16). 
We do not know which additional conditions must be imposed on X to ensure (4.9) and (4.10). 

Example 4.9. We consider a stochastic diffeomorphism as introduced by Blanc, Le Bris and Lions 
[13]; that is a field $ : R" x fl R" such that $(•, a;) is a diffeomorphism for ¥ a.e. lo Cz fl satisfying 

ess-supj^„xQ]] V3>(a;, a;)]] < M < +oo, and ess-infR^xo det V3>(a;, a;) > > 0. 
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Recalling Example 3.6 (see also [L!, Proposition let Aj(ijj) = {$(eji, Lj)}i£zn, then P a.s. in Q. 

it holds (zyA/i-"/2)ej < rA^.(w) < i?A, (w) < Msj and 



\Ju / 

By choosing Vj such that rj/sj 7 > P a.s. in Q,, we have 9 — 7" J"j^ det V3>~^(x, •) dx' P a.s. in 
fl. In conclusion, Aj are random Delone sets satisfying (4.8)-(4.10). 

Example 4.10. Consider as before a stochastic diffeomorphism, and assume in addition $ to be 
.stationary w.r.to (t±)±^z", that is for every i g Z", for a.e. x e R" and for P a.s. in fl it holds 

ri(tj)) = ^{x + ±,uj). 

LetAj{uj) ~ a;)}igz"; ™e have {i'M^~" /2)£j < r\.{uj) < R\.{uj) < Msj. Thus Kj are random 

Delone sets satisfying (4.8). For what (4.9) and (4.10) are concerned by [12, Lemmas 2.1, 2.2] and 
[13, Remark 1.9], P a.s. in il it holds 



Take also note that Aj satisifes (4.16) with Sj = 0. Hence, r^j is constant by Corollary 4-7, and 
moreover if we choose rj G {0,iyM^^"ej/2) such that rj/ej —^J>0, the T-limit of the energies in 
(4.14) is given by the functional 



for u G W''''P{U), Jf{u) = +00 otherwise in LP{U). 

Eventually, let us remark that stationary diffeomorphisms according to Blanc, Le Bris and Lions 
[13] satisfy the weaker condition 'S/'^(x,Ti_ui) = V$(a; + i,w). The sets of points generated by $ are 
not stationary according to (4.16) in general. 

Appendix A. 

We prove some elementary bounds on the singular kernels that were crucial in the computations of 
subsections 3.2 and 3.3. 

Lemma A.l. Let v > Q, then there exists a positive constant c{n, v) such that for every measurable 
set O C R" it holds 





Then by using (3.4) it follows 
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(i) if V € (0, n) and dist(z, (9) = 



-±^dx<c{n,v){C\0)f-^'-. (A.l) 



(ii) +oo) and dist(z, O) > 

1 



o 



\x — z\ 



dx <c{n,iy){dist{z,0))"~% (A.2) 



Proof. The lemma is an easy application of Cavalieri formula. 

Let us start with (i). Clearly, we may suppose £"(0) < +00 the inequality being trivial otherwise. 
Then, by setting s = (/:"(0)/cj„)i/" a direct integration yields 

\x - z\-''dx = C^{{x ^0:\x~z\< t-'/''})dt ^ V Z"'"""*''' C:\{x ^ O : \x- z\ < s]) 



s /'diam(O) \ rs /'diam(O) 

ds<iyujnj s'^-^-^ds + vCiO) j s-^'-^ds 
= ^^tjj:/"(/:"(0))i-'^/" + £"(0) ( -(diam(O))-'' + ( £!M^ ' ) < c(n, iy)(/:"(0))i-''/". 

n-V y \ UJn J J 

Inequality (A.2) easily follows from a direct integration. More precisely, we have 

\x-z\-^dx^^^{dist{z,0)r~\ 

R"\Bd.=t(..o)W 

□ 
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